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Abstract 



In a series of papers Amati, Ciafaloni and Veneziano and 't Hooft conjec- 
tured that black holes occur in the collision of two light particles at planckian 
energies. In this paper we discuss a possible scenario for such a process by 
using the Chandrasekhar-Ferrari-Xanthopoulos duality between the Kerr black 
hole solution and colliding plane gravitational waves. We clarify issues arising 
in the definition of transition amplitude from a quantum state containing only 
usual matter without black holes to a state containing black holes. Collision 
of two plane gravitational waves producing a space-time region which is locally 
isometric to an interior of black hole solution is considered. The phase of the 
transition amplitude from plane waves to white and black hole is calculated by 
using the Fabbrichesi, Pettorino, Veneziano and Vilkovisky approach. An alter- 
native extension beyond the horizon in which the space-time again splits into two 
separating gravitational waves is also discussed. Such a process is interpreted as 
the scattering of plane gravitational waves through creation of virtual black and 
white holes. 
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1 Introduction 



In recent years there has been important progress in our understanding of Planckian- 
energy scattering in quantum field theory and in string theory [Q-Q. Investigation 
of such scattering could provide us with better understanding of many fundamental 
issues such as creation of singularities and black-holes, loss of information, and string 
modification of theory of gravity. 

Amati, Ciafaloni and Veneziano and 't Hooft @ made a conjecture that black 
holes will occur in the collision of two light particles at Planck energies with small 
impact parameters. In this paper we discuss a possible scenario for such a process by 
using colliding plane gravitational waves. It was argued J]], |2j that at extremely high 
energies interactions due to gravitational waves will dominate all other quantum field 
theoretic interactions. One can then visualize the following mechanism of creation of 
singularities and/or black holes in the process of collision of such high energy particles. 
Each of the two ultra-relativistic particle generates a plane gravitational wave. Then 
these plane gravitational waves collide and produce a singularity or a black hole. It 
is difficult to perform calculations for such a process in a realistic situation. We shall 
discuss here an idealized picture. Any gravitational wave far away from sources can be 
considered as a plane wave. We assume that plane waves already have been produced 
by ultrarelativistic particles and then consider analytically the process of black hole 
formation under interaction of these plane waves. 

One has two dimensionless parameters in a scattering process at Planck energies fl], 
|2|. If E is the energy in the center of mass frame, then one defines the Planckian 
energy regime to be GE 2 /h > 1, where G is the Newton constant. This means that 
one can treat the process semiclassically. The other dimensionless parameter is GE/b, 
where b is the impact parameter. If one takes GE/b <C 1 , one can use the eikonal 
approximation. The elastic scattering amplitude A(s, t) in the eikonal approximation 
was found in HI- Fabbrichesi, Pettorino, Veneziano and Vilkovisky (FPW) gave 
the representation 

A(s,t) ~sjd 2 b e iqb e iIcl{s ' b) . (1.1) 

Here s and t are the Mandelstam variables, q 2 = —t. I c i(s,b) was taken to be the 
value of the boundary term for the gravitational action calculated on the sum of two 
Aichelburg-Sexl shock waves, 

I d (s,b) = Gs\ogb 2 . (1.2) 

This corresponds to a linearization of the problem and one cannot see in such an 
approximation black hole formation. It is tempting to suggest that one can use the 
FPVV approach J7| and the formula of type (1) even for small impact parameters to 
calculate the phase of the transition amplitude from a state containing two particles 
to a state containing a black hole. In such a case I c i(s,b) could be the value of the 
gravitational action for a solution of Einstein equation corresponding to an interior of 
the Kerr black hole created in the collision of plane gravitational waves. The parameters 



s and b in this case could be expressed in terms of the mass and angular momentum 
of the Kerr black hole. 

Classical collision of plane gravitational waves has been the subject of numerous 



investigations, see for example || BL [Hi, 11], and it has a remarkably rich structure. 



Here we are going to use the Chandrasekhar-Ferrari-Xanthopoulos duality between 
colliding plane gravitational waves and the Kerr black hole solution. 

The paper is organized as follows. Sect. 2 contains a discussion concerning quantum 
mechanical description of black hole creation in the scattering processes of particles. 
In Sect. 3 we discuss the boundary term in the gravitational action in the context of 
semiclassical approximation and present the boundary term in the form suitable to 
calculate non-linear effects in quantum collision of gravitational waves. Sect. 4 contains 
a necessary information about classical solutions describing colliding plane waves. In 
Sect. 5 we present the semiclassical transition amplitude for the process of creation 
of black hole in the collision of two plane waves. We show also that the transition 
amplitude for elastic scattering of two plane waves is equal to 1 in the leading order of 
the semiclassical approximation. This conclusion is based on the investigation of the 
behaviour of geodesies for the metric corresponding to elastic scattering of two plane 
waves. Its behaviour is described in Appendix. 



2 Black Hole Transition Amplitude 

We discuss here some issues concerning quantum mechanical description of black hole 
creation in the scattering processes of particles. Let us note first that black holes cannot 
be incorporated into the theory if we consider quantum field theory in Minkowski space- 
time. In fact it is obvious from Einstein equation 

R P u - \g^R = MGT^ (2.1) 

that if one has matter, i.e. a nonvanishing then one has a nontrivial gravitational 
field. This means that one has to start with a region of space-time which is flat only in 
some approximation. Then, in the process of collision one gets a strong gravitational 
field including perhaps black holes and/or singularities. Let us clarify the meaning of 
the transition amplitude from a state describing particles to a state containing black 
holes, 

< Black holes | Particles in almost Minkowski space-time > . 

By the analogy with the definition of transition amplitude in quantum mechanics 
this transition amplitude can be characterized by values of metric and others fields in 
two given moments of time (or by data on two given Cauchy surfaces, say E' and E"). 

Our starting point is the quantum-mechanical Feynman transition amplitude be- 
tween definite configurations of the three-metric h'^ and field $' on an initial spacelike 
surface E' and a configuration h"j and $" on a final surface E". This is 



< h", 0", E"\ti, 0', E' >= J e* s[9 >* ] V$Vg, 



(2.2) 



where the integral is over all four-geometries and field configurations which match given 
values on two spacelike surfaces, i.e. 

$\v = 0', 9k' = h' (2.3) 

= 4>", 9\v> = h" (2.4) 

In this paper we shall consider the transition amplitude in the semiclassical approxi- 
mation and we don't discuss here the introduction of ghosts and the definition of the 
measure Dg. 

We are interested in the process of black hole creation. Therefore we specify the 
initial configuration hi and 4>' on E' as configuration of gravitational and matter fields 
in Minkowski spacetime and we specify the final configuration h" and <fi" on E" as 
describing black hole. So, E' is a partial Cauchy surface with asymptotically simple 
past in a strongly asymptotically predictable space-time and E" is a partial Cauchy 
surface containing black hole(s), i.e. E" — J~(T + ) is non empty. To explain these let 
us recall some necessary notions from the theory of black holes [|T3|, [TJj] . 

Black holes are conventionally defined in asymptotically flat space-times by the 
existence of an event horizon H. The horizon H is the boundary J _ (X + ) of the causal 
past J~(X + ) of future null infinity X + , i.e. it is the boundary of the set of events in 
space-time from which one can escape to infinity in the future direction. 

The black hole region B is 

B = M- J~(l + ) 

and the event horizon 

H = j-(T + ). 

Consider a space that is asymptotically flat in the sense of being weakly asymp- 
totically simple and empty, that is , near future and past null infinities it has a con- 
formal structure like that of Minkowski space-time. One assumes that space-time is 
future asymptotically predictable, i.e. there is a surface S in spacetime that serves as 
a Cauchy surface for a region extending to future null infinity. This means that there 
are no "naked singularities" (a singularity that can be seen from infinity) to the future 
of the surface S. This gives a formulation of Penrose's cosmic censorship conjecture. 

A space-time (M, g^ u ) is asymptotically simple if there exists a smooth manifold M 
with metric g^ u , boundary T, and a scalar function Q regular everywhere on M such 
that 

(i) M — X is conformal to M with g^ v = Vfig^, 

(ii) tt > in M - T and O = on I with V M fi / on I, 

(iii) Every null geodesic on M contains, if maximally extended, two end points on 

J. 

If M satisfies the Einstein vacuum equations near X then X is null. X consists of two 
disjoint pieces X + (future null infinity) and X~ (past null infinity) each topologically 
RkS 2 . 

A space-time M is weakly asymptotically simple if there exists an asymptotically 
simple M with corresponding M such that for some open subset K of M including 



X, the region Mo U K is isometric to an open subset of M. This allows M to have more 
infinities than just X. 

The domain of dependence _D + (E) (.D _ (E)) of a set E is defined as the set of all 
points peM such that every past (future) inextendible non-spacelike curve through 
p intersects E. A space like hypersurface which no non-spacelike curve intersects more 
than once is called a partial Cauchy surface. Define -D(S) = L> + (E) \JD~(E). A partial 
Cauchy surface E is said to be a global Cauchy surface if D(Y0) = M. 

Let E be a partial Cauchy surface in a weakly asymptotically simple and empty 
space-time (M,g). The space-time (M,g) is (future) asymptotically predictable from 
E if X+ is contained in the closure of £> + (S) in M . If, also, J + (S) n J~(X + ,M) 
is contained in Z) + (E) then the space-time (M,g) is called strongly asymptotically 
predictable from E. In such a space there exist a family E(r), < r < oo, of spacelike 
surfaces homeomorphic to E which cover D + (H) — E and intersects X + . One could 
regard them as surfaces of constant time. A black hole on the surface S(r) is a 
connected component of the set 

B(t) = E(r) - J-(X + ,M). 

One is interested primarily in black holes which form from an initially non-singular 
state. Such a state can be described by using the partial Cauchy surface E which has 
an asymptotically simple past, i.e. the causal past J~(E) is isometric to the region 
J~(I) of some asymptotically simple and empty space-time with a Cauchy surface X. 
Then E has the topology R 3 . 

In the case considered one has space-time (M, g^ v ) which is weakly asymptotically 
simple and empty and strongly asymptotically predictable . 

E' is a partial Cauchy surface with asymptotically simple past, E' ~ R 3 . 

E" = S(t") contains a black hole, i.e. E"— J~(X + ,M) is nonempty. 

In particular, if one has the condition E' n J~ (X) is homeomorphic to R 3 (an open 
set with compact closure) then E" also satisfies this condition. 

Strictly speaking one cannot apply the standard definition of black holes to the 
case of plane gravitational waves and we need a generalization. One defines a black hole 
in terms of the event horizon, J~(X + ). However this definition depends on the whole 
future behaviour of the metric. There is a different sort of horizon which depends only 
on the properties of space-time on the surface E(r) [|13|. Any point in the black hole 
region bounded by r = 2m in the Kruskal diagram represents a trapped surface (which is 
a two-dimensional sphere in space-time) in that both the outgoing and ingoing families 
of null geodesies emitted from this point converge and hence no light ray comes out 
of this region. A generalization of the definition of black holes in terms of trapped 
horizon has been considered in [?, ^TJ. A generalization of the standard definition of 
black holes to the case of nonvanishing cosmological constant was considered in [?]. 

We discussed the transition amplitude (propagator) between definite configurations 
of fields, < h", (f>", E"|/i', </>', E' >. The transition amplitude from a state described by 
the wave function ^/ in [h',(p'] to a state ^/ out [h", 4>"] reads 



< \$/ out m in >— 



(2.5) 



f y out [ti',(f)"} < h",<j>", £'>', 0', £' > ¥ n [ti,(f>')VtiV<f>'Vh"V<f>". 

One can take for example the state ^ m = ^ l p lP2 [h', 0'] as a Gaussian distribution 
describing a state of particles with momenta p% and P2 and take ^> out as a wave function 
describing a state of black hole. Recently Barvinski, Frolov and Zelnikov have suggested 
an expression for the wave function of the ground state of a black hole 

3 Boundary term in Gravitational Action and Semi- 
classical Expansion 

In this section we discuss an approximation scheme forcalculating the transition ampli- 
tude following to the FPVV approach J7|. The gravitational action with the boundary 



term has the form [15| 



= TpTt / R V=9d 4 x - — / KVhd 3 x (3.1) 
167rG Jv 8nG Jav 

Here V is a domain in space-time with the space-like boundary dV, h is the first 
fundamental form and K is the trace of the second fundamental form of dV . 



The case of null surfaces was considered in \T&. We shall write a representation of 



the boundary term suitable for quantum consideration. The action is 

S[g] = -j^ (J y d 4 x^-g R(g) + jf d^x^g V„/"G/)) , (3.2) 

where 

F(g) = g aP g" v d u g aP - g^g Pv d u g aP . (3.3) 
Supposing that the boundary is described by equation 

a(x) = (3.4) 

one gets the action in the form 

%d] = j y dV=0 R{g) + J y d'x^g~ l 5(a(x)) /"V M <r). (3.5) 

The linearization of the action ( |3.2|) leads to the action in the FPVV form [|7|. 

Let us show that the presence of the boundary term in the action ( |3.2j ) is necessary 
for the selfconsistency of the semiclassical expansion. To perform semiclassical expan- 
sion one expands the metric g around a classical solution g c i of the Einstein equation 
so that g = g c i + Sg, 

S[g] = S[g d + 5g] = S[g cl ] + S'[g cl ]5g + ±S"[g cl ](5g) 2 + (3.6) 
where g c i matches the given Cauchy data on surfaces £' and i.e. 

9d Is' = hf, g d | s » = h" (3.7) 



In this case 5g\x' = <5g|s» = 0, but we cannot guarantee that VSg\s' = V<5g|s" = 0. To 
ensure that terms linear on V5g drop out from expression (|3.6| ) (otherwise we cannot 
perform semiclassical expansion) one has to integrate by parts. One has 

%d + Sg] = (| d 4 xV=9li R(9d) + J d 4 xV=9li V^GfcOH (3.8) 

+ / d i x^g~ l (g^V 2 5g, u - V»V»5g, u ) + V^PSg^gT ~ ^<W)) 

+ second order terms. 
The linear terms coming from the Hilbert-Einstein action can be put in the form 

sTc^VaVpSg^ - = V,(g a ^5g aP - g^Sg^) (3.9) 

Notice that on the RHS of the last relation covariant derivatives can be replaced by 
partial derivatives because on the boundary 5g = 0. Therefore one finds that terms 
linear in 5dg coming from Hilbert-Einstein action cancel similar terms coming from 



full divergence and the action ( |2.2| ) admits the expansion ( |3.6| ). 

Taking into account that the value of the Hilbert-Einstein action on the classical 
solution is equal to zero one finds that the full action for a solution of Einstein equation 
is reduced to the second term in (|3.8|) which can be reduced to boundary term. 

The transition amplitude ( |2.2| ) in semiclassical approximation is 



< h", E"\ti, £' >= Zexp % -S d (3.10) 

n 



where 



%d] = -ttt^ / d'x^g~ d 6(a(x)) f»V »<j, (3.11) 
lo7rG J 

f is given by equation (FO) and 



Z=( , 1, J 1/2 



det S"(g 



cl 



We assume here that there is only one solution of classical equation of motion with 
given boundary conditions. 

4 Colliding Plane Gravitational Waves 
4.1 Szekeres line element 

In this section we present a solution of Einstein equations describing colliding plane 
gravitational waves. There exists a well known class of plane-fronted gravitational 
waves with the metric 



ds 2 = 2dudv + h(u, X, Y)du 2 - dX 2 - dY 2 



(4.1) 



where u and v are null coordinates . In particular the gravitational field of a particle 
moving with the speed of light is given by the Axelburg-Sexl solution. The metric has 
the form 

ds 2 = 2dudv + E log(X 2 + Y 2 )5(u)du 2 - dX 2 - dY 2 (4.2) 

and describes a shock wave. It is difficult to find a solution which describes two sources. 
An approximate solution of equation (j2.1|) for two particles as the sum of solutions each 
of which describes one particle was considered by FPVVjTj. This approximation de- 
scribes well the scattering amplitude for large impact parameter, but does not describe 
non-linear interaction of shock waves which is dominate in the region of small impact 
parameter. To analyze non-linear effects we will, instead of dealing with shock wave, 
take a simple solution of Einstein equation, namely we will take plane gravitational 
waves. In some respects one can consider plane wave as an approximation to more com- 
plicated gravitational waves, in particular shock waves. This solution in some sense can 
be interpreted as an approximation for a solution of Einstein equation in the presence 
of two moving particles. Collision of two ultrarelativistic black holes was considered 
by D'Eath JlTj. If our particles are gravitons then there are not sources corresponding 



to matter fields in Einstein equations. Note also that plane gravitational waves are 
produced by domain walls, see |3Q. 

A particular class of plane waves is defined to be plane-fronted waves in which 
the field components are the same at every point of the wave surface. This condition 
requires that h(u,X, Y) is a function with a quadratic dependence on X and Y. One 
can then remove the dependence of h on X and Y altogether by a coordinate change. 
Solutions of Einstein equations describing collisions of plane gravitational waves were 
first obtained by Szekeres and Khan and Penrose |§. Chandrasekhar, Ferrari and 
Xanthopoulos || |10[ have developed a powerful method for obtaining such solutions 
by using a remarkable analogy ("duality") with stationary axisymmetric case which can 
be reduced to the investigation of the Ernst equation. For a review see the Griffiths 



book 11 



We will use the coordinates (u,v,x,y). We assume that throughout space-time 
there exists a pair of commuting space-like Killing vectors £1 = d x , £2 = d y . The 
Szekeres line element has the form 

ds 2 = 2e~ N dudv - e' u (e v cosh Wdx 2 - -2 sinh Wdxdy + e v cosh Wdy 2 ), (4.3) 

Here N, U, V and W are functions of u and v only. 

We illustrate in Fig.l the two-dimensional geometry of plane waves. Space-time 
is divided into four regions. The region I is the flat background before arriving plane 
waves. The null hypersurfaces u = 0, and v = are the past wave fronts of the incoming 
plane waves 1 and 2. The metric in region I is Minkowski. Regions II and III represent 
incoming plane waves which interact in region IV. Colliding plane gravitational waves 
can produce singularities or Cauchy horizons in the interaction region |], [1^. [20|. [21 



A solution is undetermined to the future across a Cauchy horizon || 0, [21], P3| . We 
shall discuss two simplest extensions. 

In particular one can get an interiour of the Schwarzschild solution in the interac- 
tiong region IV. There are two types of colliding plane waves solutions corresponding 



to the Schwarzschild metric. The first one creates the interior of the blach hole with 
the usual curvature singularity. In this case incoming plane waves have curvature sin- 
gularities already before collision. In the context of the Planck energy scattering it 
seems more natural don't have curvature singularities already for free plane gravita- 
tional waves. Therefore we mainly will be discussing another type of solutions when one 
getsin the interaction region the interior of the Schwarzschild white hole. The maximal 
analytic extension of this solution across its Killing-Couchy horizon leads to creation 
of a covering space of the Schwarzschild black hole out of collision between two plane 
gravitational waves. An alternative interpretation of this solution is the creation of 
the usual Schwarzschild black hole out of the collizion between two plane gravitational 
waves propagatiing in a cylindrical universe. There exist also a time-reversed extension 
21 ] including covering space of the Schwarzschild exterior and part of black hole, and 



giving two receding plane waves with flat space between. We will interprete this as 
scattering of plane waves on the virtual black hole. 

4.2 Ernst's equation for colliding plane gravitational waves 

Vacuum Einstein equations 



- = (4.4) 



for the metric (4.3) have the form 



U uv = U U U V (4.5) 

2U VV = U 2 + W 2 + V 2 cosh 2 W - 2U V N V (4.6) 

2U UU = U 2 + W 2 + V 2 cosh 2 W-2U u N u (4.7) 

2U UV = U u V v + U v V u -2(V u W v + V v W u )tanhW (4.8) 

2N UV = -U u V v + W u W u + V u V v cosh 2 W (4.9) 

2W UV = U u W v + U v W u + 2V u V v smhWcoshW (4.10) 



One can reduce the system of equations ( |4.5|) -( [O"0|) to a solution of the Ernst equation 
for the complex valued function 

Z = X + i\ (4.11) 

where 

X = e~ v sinh W, X = e~ v tanh W. (4. 12) 



The line element becomes (JO) 



ds 2 = Ae~ N dudv - e~ u [xdy 2 + -{dx - Xdy) 2 ], (4.13) 
First note that equation (|4.5|) can be integrated to give 



e- u = f(u) + g(v), 



(4.14) 



where f{u) and g(v) are arbitrary functions. We fix a gauge by using / = f(u) and 
g = g{y) as new coordinates instead of u and v and then we change variables from / 
and g to \i and 77 such that 

/ + 2=(1-^)V 2 (1-?7 2 ) 1/2 , f-9 = -M- (4-15) 
One has the Ernst equation 

(Z + Z)[((l - /.X),, - ((1 - ^Z,,),,] = 2[(1 - /i 2 ^ - (1 - r?)Zl\. (4.16) 

If a solution Z of the Ernst equation ( |4.16|) is found then the metric functions V and 
W are given by 

e v = (ZZ)-\ sinhiy = -i Z - Z . (4.17) 

Z -\~ z 

U is given by e u = (1 — / u 2 ) 1 / 2 (l — 77 2 ) 1 / 2 and iV is defined by the relation e~ N = 
f'd'if + g)~ 1 ^ 2 e~ B . To find N one has to integrate the following equations 

Bj = -2(f + g)(l- EE)- 2 EjE f , B, g = -2(f + g)(l- EE)- 2 E, g E, g , (4.18) 

where E is the Ernst function related with Z, E = [Z — 1)(Z + 1) _1 . 



4.3 Chandrasekhar-Ferrari-Xanthopoulos duality 

The Ernst equation ( J4.16 ) is invariant under the Neugebauer-Kramer involution. If 



Z = x + *A is a solution of Q4.16|) then 

Z = $ + W 



is also a solution of (|4.1( 



where $ and \l/ are given by equations 
* = (1 



^^(l-^V^-i. 



— 7F~ A .»? 



(4.19) 



(4.20) 
(4.21) 



Notice here that for stationary axisymmetric space-time one has the same Ernst equa- 
tion (|4.16|) , whose solution Z = \& + 2$ is considered as a potential for the metric 
functions x an d A that are obtained from ( |4.21| ). Take 



E = p/i + iqrj. 
This is a solution of Ernst's equation if 



p 2 + q 2 = 1 



(4.22) 



(4.23) 



If one considers ( [4.22[ ) as containing the metric functions, E = x + *A, it leads to the 
Nutki-Halil solution describing colliding gravitational waves with non-aligned polar- 
ization. However, when the same function (|4.22|) is considered as a potential, it then 



leads to the Kerr solution. Therefore a region of the Kerr space-time can be also con- 
sidered as a solution describing plane waves. This remarkable observation we call the 
Chandrasekhar-Ferrari-Xanthopoulos duality. The Ernst potential ( [4.22 ) gives 

Z = * + i$ = i±^±M (4.24) 
1 — pfi — iqrj 

and one gets the following solution of equations ( |4.20|) and ( |4.21| ) 

X = (l-^ 2 ) 1/2 (l-^ 2 ) 1/2 f (4-25) 



2q 1 (1-7^(1-^ 
p 1 + p 1 — p 2 yU ? — <7 2 ?7 2 

Here 

X = (1 - P/ u) 2 + gV , y = 1 - pV - gV (4.27) 
For the Ernest potential the metric has the form 

ds 2 = X(-^- 2 - -^L) - (1 - p 2 f' 2 {l - V 2 )^\xdy 2 + ~{dx ~ \dyf\. (4.28) 

1 — if 1 — p 2 X 

Now let us take 

p=-- — , q = — 1 m>a (4.29) 

m m 

and introduce the new coordinates (t, r, 9, 0) instead of (p, 77, x, y) by putting 

r — m 

^=7"^ 2TT72"' V = cos6 4.30 

r- 2q \plm 

t = -V2m(x- ^j—j), t = {m2 _ a2)l/2 y (4.31) 

Then the metric ( |4.28| ) will take the form of the Kerr solution 

2m 2 ds 2 = (l-^^)dt 2 sin 2 9dtd(j)-(r 2 +a 2 - — ™' r ) sin 2 6d(j) 2 -p 2 (-^dr 2 +d6 2 ) ) 
p 2 p 2 p 2 A 

(4.32) 

where 

p 2 = r 2 + a 2 cos 2 9, A = r 2 - 2mr + a 2 (4.33) 

The coordinates must satisfy the inequality \rj\ < p < 1. This implies that — (m 2 — 
a 2 ) sin 2 6 < A < 0, which means the region of the Kerr spacetime that is inside the 
ergosphere. 

To describe the colliding plane gravitational waves producing the interior of the 
ergosphere in the Kerr spacetime it is convenient to rewrite the metric (|4.28|) in terms 
of the (u f , v') coordinates related with (rj, p) by the following relations 

rj — sin(w' — v'), p = sin(w' + v'), (4.34) 



These (u',v') are related with (u, v) in (|4.3|) by u = sin-u', v = sin?/. For simplicity of 
notations we will omit ' in ( 4.34j) . We have 

ds 2 = AX(u, v)dudv — Q(u, v)[x{u, v)dy 2 H - -(dx — X(u, v)dy) 2 ], (4.35) 

where 

X(u, v) = (1 — psin(u + v)) 2 + q 2 sin 2 (u — v), Q(u, v) = cos(m + v) cos(m — v), (4.36) 
2q 1 — sin(u + v) 



X(u, v) 



-((p + 1) sin (u — v) + psin(u + v) — 1) 



(4.37) 
(4.38) 



l+p Y(u, v) 

Y(u, v ) = 1 — p 2 sin 2 (w + v ) — q 2 sin 2 (w — v). 

In (|Q5D 

0<m<tt/2, 0<v<ir/2, v + u<tt/2 (4.39) 
To extend the metric (|4.35|) outside of region ( |4.39|) one uses the Penrose-Khan trick 

u -> u9(u), v -> v6(v), (4.40) 

« (*) = I n' ! ^ (4-41) 



and substitutes in ( |4.35| ) 



Fig.l illustrates this metric. The region I (u < 0, u < 0) is Minkowskian. Regions II 
and III contain the approaching plane waves with the following metrics 



where 



(ds 11 ) 2 = AX(u)dudv - n(u)Mu)dy 2 + —-(dx - \(u)dy) 2 ] 

xw 



(ds 111 ) 2 = AX{v)dudv - tt(v)lx(v)dy 2 + — -(dx - \{v)dy) 2 ] 

X\P ) 

X(u) = (1 — 2p sin u) + sin 2 u, Q(u) = cos 2 u, 

1 4 



X(u) = 2g(sinw 



l+p 

u < 7r/2, u < 0, 



In (|442D 
and in flOj) 

it < 0, u < 7r/2, 
The region IV is the interaction region with metric ( f4.35|) . 



(4.42) 

(4.43) 

(4.44) 
(4.45) 

(4.46) 
(4.47) 



4.4 Plane Gravitational Waves for Zero Impact parameter 



In this subsection we discuss a particular example of two colliding plane waves, namely 
p = — 1, ||21|| . The choice is motivated by the fact that the colliding plane gravitational 
waves produce in the interaction region a space-time that is isometric to an interior of 
the Schwarzchild solution. The metric is given by 

ds 2 = 4m 2 [l + sm(u9(u)) + v9{y)]dudv 

- [1 - sm{u9(u)) + v9(v))[l + sm(u9(u)) + v9{v))~ 1 dx 2 (4.48) 

-[1 + sin(u0(u)) + v9(v)} 2 cos 2 (u9(u)) - v9(v))dy 2 , 

where u < tt/2, v < tt/2, v + u < ir/2. 

Fig.l illustrates this solution of the vacuum Einstein equations. The background 
region I describes a region of space-time before the arrival of gravitational waves and it 
is Minkowskian. Two planes waves propagate from opposite directions along the z-axis. 
Regions II and III contain the approaching plane waves. In the region IV the metric 
( [4.48 ) is isomorphic to the Schwarzchild metric. To see this one can make the following 



change of variables from "plane waves" coordinates to Schwarzchild coordinates, 

(u,v,x,y)->(t,r,9,<l>) (4.49) 

defined by, 

r — m[l + sin(w + v)), t = x, 9 = n/2 + u — v, (f) = y/m, (4.50) 
or to Kruskal coordinates t,(,9,4> 

t = —a(r) cosht/4m, ( = —a{r) sinht/4m, (4-51) 

where 

a(r) = (1 - r/2m) 1/2 e r/4m . (4.52) 

Then one gets 

ds 2 = 6 — e - r l 2n \dT 2 - dC) - r 2 {d9 2 + sin 2 9d<\> 2 ) 
r 

The section of the region IV bounded by x = 0, y = corresponds to segment in the 
Kruskal diagram and the section of the region IV by the plane x = xo, yo = corre- 
sponds to the hatched region in the Kruskal diagram Fig.5. The lines corresponding to 
r = 2m (horizon) apart from the point (r = 0, £ = 0) correspond to the infinite value 
of x-plane wave coordinate. 

The above metric in the (u, v) plane can be extended beyond the event horizon 
u + v = tt/2 in one of two ways. 

The first possibility, shown in Fig. 6 consists in reflecting along the line u + v = it/ 2. 

The second one involves in gluing to the horizon the whole upper-half part of the 
Kruskal diagram. 



Both extensions are solutions of Einstein equations. From the general discussion 
in Section 3 there is a-priori non-zero probability to get a finite state corresponding 
to a black hole or two outgoing plane waves. In the next section we will calculate the 
probabilities for these processes in the semi-classical approximation. As was mentioned 
in the Sections 2 and 3 to define correctly the finite Cauchy surface we have to know 
how the null geodesies look like in space-time with the metric ( |4.48| ) extended as shown 
in Fig. 6. 

For the usual black hole solution, i.e. the regions IV and IV are extended up 
to the boundary r = and there are two asymptotically free regions, namely A and 
A in fig. 5. In this case some geodesies cross the lines r = ±2m and may not return 
to region VI . So, it is important to know if it is possible for geodesies starting from 
Minkowski region I to end up in the region A or A . A proof of this fact is given in 
Appendix. 



5 Semiclassical Transition Amplitude 

In this section we study the transition amplitude 

< BH\2pw > (5.1) 

from a state \2pw > of two plane gravitational waves to a state < BH\ containing 
black-hole, and the transition amplitude 

< 2pw\2pw > (5.2) 

from two plane gravitational waves back to two plane gravitational waves. We can 
consider these two amplitudes as amplitudes of two independent channels. 

To find these transitions amplitudes in the semiclassical approximation according 
to ( |3.10| ) we have to evaluete the boundary term on the classical solution interpolated 
between two plane gravitational waves in the initial state and black hole in the final 
state, and on the classical solution interpolated between two plane gravitational waves 
in the initial state and two plane gravitational waves in the final state. 

Let us start from consideration the transition amplitude 

<2pw,WH\2pw> (5.3) 

from two plane gravitational waves to two plane gravitational waves and white hole. 
The corresponding classical solution is shown on fig. 3. In this case the initial Cauchy 
surface £' crosses regions I, II and III and the final surface crosses regions II, IV and 
III (Fig. 7). Suppose that the equation ( |3.4| ) for the surfaces E' and E" has the form 

£' : a = u - v - w' = 0, w' Q < (5.4) 

E" : o~ = u — v — Wq = 0, < w'l < tt/2 (5.5) 



In this case the other parts of the boundary defining the boundary term are given by 
equations 

B 2 : u - vr/2 = 0, w' -7I-/2 < v < w'q - tt/2; B' 2 : u = 0, w' - tt/2 < v < (5.6) 

fi 3 : u - tt/2 = 0, w' - tt/2 < u < w% - tt/2; B' 3 : w = 0, w£ - tt/2 < u < (5.7) 

To calculate the phase factor let us give expressions for / M on the boundary. In 
Minkowski space we have / M = 0. In all regions f x = f y = 0, moreover 

r\n = o, rim = o (5.8) 

So the value of the action (|3.11|) for the two plane waves solution ( f4.42|) , (^4.43|) is 
reduced to the sum of two terms each of which represents a contribution from the 
Cauchy surface £' and respectively, 

S(g% pw) ) = + (5.9) 

where 

5 s ' = S 2 + S' 3 , 5 s " = S% + S' 3 f + S'l, (5.10) 

see fig. 7. So 

< 2pw\2pw, WH >= Zexp{-—^—(S' 2 + S' 3 - S' 2 ' - S» - S'l)} (5.11) 

here Z is the one-loop contribution. 

Taking into account the special form of the metric one can write f w in regions II, 
III and IV respectively in the following form 

y/=gr\n = ^(u)d u \n[X(u)n 2 (u)} (5.12) 

V=jjr\iu = V(v)d v \n[X(v)n 2 (v)} (5.13) 

y/=gr\iv = z)d™ HX(w, z)n 2 (w, z)] (5.14) 

W — U + V , z = v — u. 

The representation (]5.12|) was obtained in the following way. Metric ( [4.35 ) in (w, z, x, y) 
coordinates admits the representation 

with the diagonal matrix g a p, 

9ap = ( 9 q W ® ^ •, 9ww = ~9zz = X (5.16) 

and Qij of the form 

0<j=%y (5.17) 



where 

/ i A 



being an element of the group SL(2). From equation ( |3.3| ) it follows 

f w = 9 ww (g zz d w g ww + g ij d w9ij ) (5.19) 
Taking into account that the trace trg~ 1 dg = for the group SL(2), we get 

f w = g ww {g zz d w g ww + 2tt- l d w tt), (5.20) 



from which follows the representation (|5.14 ). 
Simple calculations give that 

1 rL/2 rL/2 rO rw' 

S' = / dx dy / dv5{u + v-w')tt(u)d u \n[X{u)tt 2 {u)} 

167TG J-L/2 J-L/2 Jn/2 Jw' -n/2 



1 f° 

— — -L 2 / du cos 2 uO u [ln((l — p sin u) 2 + q 2 sin 2 u) cos 4 u] 
16n& Jtt/2 



(5.21) 



1 ; 2L'(/V 1 -»'X''-?> -4 /•«*.) 



16vrG Vi 1 - 2up + u 2 



-3 — 2p + 2g In 2(1 — p) + Apq arctan ) 

167rG 1 — p 

where L is the cut-off in the x — y-plane, 

S >> = S> 2 + SI (5.22) 

where 

L 2 r° 

So = — / du cos 2 ud u [\n((l — psinu) 2 + q 2 sin 2 u) cos 4 w] = (5.23) 
lo7rG Jw% 

X q sin uu 

■ [—3 sin 2 w'q — 2p sin w'q + 2q 2 ln(sin 2 w' Q ' — 2p sin w' ' + 1 ) + Apq arctan ■ 



L , v „ „ , , -j-^ „ . j 

167rG l-psmwo 

The part of the boundary which cross the white hole region gives the following 
contribution 

1 rL/2 rL/2 rw' 

S'l = / dx dy dzn(w,z)d w ln\X(w,z)n 2 (w,z)}= (5.24) 

167rG J-L/2 J-L/2 J-w' 

L 2 r~ w o , . — 2p( 1 — p sin tfn ) cos 2 w'J, cos z ,. , 

' d«[— ^ 5 °\ „ ,° 2 2 sin < cos z] 



167rG Jw'j (1 — psin w'q) 2 + g 2 sin z 

In particular, for Wq = 7r/2 only the second term in the RHS of ( |5.24|) survives and we 
get 

_™ .... L 2 



Therefore the amplitude describing the transition from two plane waves to wormhole 
in the semiclassical approximation is given 



< WH\2pw >= Zexp{Us' 2 + S' 3 + S'{)} = (5.25) 

Zexp{ — — — — [— 1 — 2p + 2q 2 In 2(1 — p) + Apq arctan ]} 

StyCtTl 1 — p 

To get the semiclassical answer for the transition amplitude < 2pw\2pw > we have 
to consider the classical solution which describes the metric prolongated as shown in 
Fig. 6. We see from equations ( |5.21|) that S' 2 does not depend on w' and therefore for 



the transition amplitude describing the elastic scattering of two plane waves the phase 
factors is zero, and 

< 2pw\2pw >= Z (5.26) 



CONCLUSION 



We have considered a possible mechanism for black hole creation in the collision of 
two light particles at planckian energies. Many questions deserves a further study. 
In particular there is an important question how to relate the momenta of colliding 
particles with characteristics of colliding plane waves. It is known that a plane wave 
doesn't lead to polarization of vacuum p6[. It was shown that a plane wave is an 



exact string background |27L gg, |29], |30[ and there is a duality with extreme black hole 



TL\. One can conjecture that a dilaton gravity analogue of the CFX duality between 
colliding plane waves and non-extremal black holes discussed in this paper could lead 
to a corresponding string duality. 
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A APPENDIX 



Geodesies for colliding plane gravitational waves 
with zero impact parameter 

Investigation of geodesies in the metric describing two colliding plane waves has been 
considerad in Ref. []. As has been mentioned above we are interested in behavior of the 
null geodesies starting from the surface £' (fig. 7). Here we prewsent the explicite form 
of geodesies in the metric ((O^) describing two colliding plane waves for zero impact 
parameter. The case of metric Q4.35p is more complicate but there is a raison to believe 
that the behavior of geodesies for the case of zero impact paremeter reproduces the 
essential feature of the general case. 

Suppose that a geodesic starts in the Minkowski region I with initial momenta p M . 
We investigate its trajectory upto the region I and we express the final 4-momentum 
in terms of its initial values in region I (p x ,p y , Pu,Pv)- The geodesic starting in the 
region I at the point (uq, Vq, Xq, jJq) crosses the boundary between the regions I and II 
at the point (0, v^ 1 , Xq 1 , y^ 1 ) determined by the following equations 

tf-< = ~< (A.l) 

Pv 

II I _ ^Px m2 I /. n\ 

Pv 

vi 1 - vl = 2j ^< (A.3) 

Pv 

It is well known [?] that to get the explicit form of geodesies 

5? + K p x u x p = 0, (A.4) 
it is convenient to use the first order equations 

m x )1 = g^Pv, (A. 5) 

where p v = 4^ and the action satisfies the Hamilton- Jacobi equation 



= < (A.6) 

dx^dx u 



For the metric ( 4.48| ) in region II three components of p^, namely, p v ,p x ,p y , are integrals 



of motions, and the Hamilton- Jacobi action S in this region has the following form in 
plane wave coordinates 

S 11 = dup u + p v v + p x x + p y y (A. 7) 

Juo 

where u-component of momentum is given by 

m 2 r (1 + sinw) 4 „ 1 „ . „ „ 

Pu = — {- ^^Pl + — ^^ + (l + sm M ) 2 m 2 }, (A.8) 

p v cos z u cos z u y 



(mass-shell condition). 

Equation (|A.5|) in can be solved explicitly, 



2m 2 p x r u , (1 + sin it) 4 2m 2 p x , , rr . . . . 

x-Xo J = — / d«- — '— = — [8(tanM-tan^ ) + (A.9) 

p v J Uq 1 COS U p v 

XX X X 5 

8( jj ) + 4(cosm — cosMq 7 ) + -(sin2u — cos2uq 7 ) — —(it — it 7 /)], 

COS U COS U ^ 

TT 2m 2 p v f u du 2m 2 p Vr TT . .. 

2/ - % = — / — 5- = — t tan u - tan «o ] ( A - 10 ) 



v - « " = -**.{x - 4 1 ) y» + (A.ll) 

Pv Pv 

3 1 1 

mg[-(t! — it 7 /) — - costt(l + sin it) + - cosUq / (1 + sin tig 1 )], 

Here the initial data it 7 /, Vq 1 , Xq 1 , y^ 1 , P^PviPxiPy f° r the region II are given by the 
final data in the region I at the boundary between regions I and II, and therefore, 

The action S in regions IV and IV in the Schwarzschild coordinates has the fol- 
lowing well-known form 



S = -Et + l (f) + J ] jL 2 - -^- e d9 + J sjE 2 - (1 - ^)(mg + ^) T ^, (A.12) 

where I, L and E are three integrals of motions related to the final data for the region 
II by 

E = p x (A.13) 
/ = mp y (A. 14) 



m 2 p 2 y m 2 (l + sinixp) 4 2 _p, 2 

COS^ Uq p v COS z Uq COS z Uq 



y _ + { _[ V " p 2 + + m 2 (1 + sin%)2] _ p j2 (A _ 15) 



Here it determines the point (it , 0) in the (u, v) plane at which the trajectory crosses 
the boundary between II and IV. Let us discuss for simplicity the case p x = 0. In this 
case for null-geodesics the coordinates of the point Uq V , 0, x^ v , y^ v , u 1 ^ = uq are given 
in terms of the initial data in the Minkowski region by 

p 2 1 

m = arctan(- — v —Vq + — Uq) (A. 16) 

2mpy 2m 

4 V = 4 (A- 17) 

= y z + 4 + — <■ (A.i8) 

P 2 y Pv 



L 2m 



Equations for geodesies in the regions IV and IV for E = are 

i = 
m§r 

m <fi 
m 9 



- 1 



-I 



r 2 sin 2 9 
-\jL 2 - I 2 / sin 2 9 



(A.19) 
(A.20) 

(A.21) 
(A.22) 



The trajectory in the (r, 9) plane is determined by equations ( |A.20| ) and ([A. 22]) with 
the initial data 

r\ e =e =ro, r = m(l + sinw ). (A.23) 



Integrating, one finds 
cos 9 



cos 6*0 



arcsm 



— arcsm 



1 - l 2 /L 2 



1 - l 2 /L 2 



arcsin(l — r/m) — arcsin(l — r /m) (A. 24) 



Using the above solution we can find the point (tt/2,v ) in the (u, v) plane where the 
trajectory crosses the boundary between regions IV and II . Indeed, on the boundary 
between regions IV and II we have 



9> = 9\ 



«=tt/2, v=v — 7T — Vq 



and 



r n = r 



u=7T /2, v=v =m(l + cos w ). 



(A.25) 



(A.26) 



Equation ( |A.24j ) also gives the relation between r' and 9' Q . Substituting in that re- 
lation the expressions of r and r' Q in term of u and v ( equations (??) and ( |A.26| ) 
respectively) and 9 = n/2 + u and (|A.25|) we get 



COSfQ 



1 - 



sm uq 



smwo 



1 - l 2 /L 2 \ 1-P/L 2 ^ Jl-P/L 2 \ 1-P/L 



1 - 



COS z Vq 



cos(f + Mo) (A. 27) 



Let us look at the geodesies in the (r, () plane of the Kruskal coordinates. From 
equations (|A.19HA.22|) and (|4.51|) it follows that 



which yields 



where 



dC = C 
dr t ' 



, Co 
Cr = — r, 

To 



Co = -a(r(t ))sinht /4m, 



(A.28) 

(A.29) 
(A.30) 



Taking a suitable initial data in the Minkowski region I we can set Co = (this choice 
corresponds to motion in the x = plane). All geodesies then pass through the 
wormhole in the Kruskal diagram. 

In region II the geodesies are similar to those in the region II with p v ,p u inter- 
changed, i.e. 

S 11 = / dup u + p v v + p x x + p y y (A.31) 

JUq 

where the f-component of momentum is given by 

m 2 r (1 + cost> ) 4 9 (1 o ,9 . . 

Pv = ■ 2 P i + -r^-Pl + 1 + cosv) 2 m 2 }, (A.32) 

p u sin v sin v y 

p u can be found from the equation 



~2 4^ - - 2 . 4 + P» -P« + 2m ( — I r-r— ) (A.33) 

p* v cos* m j>„ sin v o cos z m sin u 

which follows from junction conditions on the boundary between II and IV and II 
and IV regions. u and v in ( |A.33| ) are defined by equations ( JA.16Q and ( |A.27| ). 

From above consideration follows that all null geodesies starting from the Cauchy 
surface £' reach the surface S" or end up on the null surfaces u = tt/2,v < 0; and 
v = tc/2,u < 0. 



(ds 1 ) 2 = Am 2 dudv — dx 2 — dy 2 , 
(ds 11 ) 2 = 4m 2 [l + smu}dudv — cos 2 m[1 + sinu]~ 2 dx 2 — cos 2 w[l + sinu] 2 dy 2 
(ds 111 ) 2 = 4m 2 [l + smv]dudv — cos 2 v [1 + sinv]~ 2 dx 2 — cos 2 t>[l + sinv] 2 dy 2 
(ds IV ) 2 = Am 2 [1 + sin(w + v)]dudv — cos 2 (w + v)[l + sin(w + v)]~ 2 dx 2 

— cos 2 (w — v)[l + sin(w + v)] 2 dy 2 

Figure 1: (u,v) plane wave coordinates 




Figure 2: (x, u, v) plane wave coordinates 



IV 




Figure 3: Kruskal coordinates 



(ds 1 ') 2 = 4m 2 dudv — dx 2 — dy 2 , 
(ds 11 ') 2 = 4m 2 [l + cosv]dudv — sin 2 -y[l + cosv]~ 2 dx 2 — sin 2 -i;[l + cosv} 2 dy 2 
(ds 111 ') 2 = Am 2 [1 + smu]dudv — sin 2 -u[l + cosu]~ 2 dx 2 — sin 2 w[l + cosu] 2 dy 2 
(ds IV ') 2 = 4m 2 [l + sin(u + v)]dudv — cos 2 (w + v)[l + sin(w + v)]~ 2 dx 2 

— cos 2 (w — v)[l + sin(w + v)] 2 dy 2 

Figure 4: Extended plane wave coordinates 




Figure 5: Initial and final Cauchy surfaces 



